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In this paper we study the polarization asymmetries of the final state lepton in the radiative dilep- 
tonic decay of B meson {Bs £'^£~'y ) in the framework of Minimal Supersymmetric Standard 
Model (MSSM) and various other unified models within the framework of MSSM e.g. mSUGRA, 
SUGRA (where condition of universality of scalar masses is relaxed) etc. Lepton polarization, in 
addition of having a longitudinal component {Pl ), can have two other components, Pt and Pjv , 
lying in and perpendicular to the decay plane, which are proportional to me and hence are significant 
for final state being /i+ fi~ or r"*" r~. We analyse the dependence of these polarization asymmetries 
on the parameters of the various models. 

PACS numbers: 13.20He,12.60.-i,13.88+e 



I. INTRODUCTION 

Flavor Changing Neutral Current (FCNC) induced B-meson rare decays provide a unique testing ground for Stan- 
dard Model (SM) improved by QCD corrections via Operator Product Expansion ( for a review and complete set of 
references see [Q). Studies of rare B decays can give precise information about various fundamental parameters of 
SM like Cabbibo-Kobayashi-Maskawa (CKM) matrix elements, leptonic decay constants etc. In addition to this, rare 
B decays can also give information about various extensions of SM like two Higgs doublet model (2HDM) |§, ||, ^, ||] 
, Minimal Supersymmetric Standard Model (MSSM) §, §, |, § |l0[ |ll], |l|| etc. After the first observation of 
the penguin induced decay B Xs^y and the corresponding exclusive decay channel B K*j by CLEO , rare 
decays have begun to play an important role in particle physics phenomenology. 

Among the rare B decays, Bg — > £^£~^ {£ = e,fj,,T) are of special interest due to their relative cleanliness and 
sensitivity to new physics. They have been extensively studied within SM 17, |l^ and beyond [^. In the mode 



Bg £~^£~^ , one can study many experimentally accessible quantities associated with final state leptons and photon 
e.g. lepton pair invariant mass spectrum, lepton pair forward backward asymmetry, photon energy distribution and 
various polarization asymmetries (like longitudinal, transverse and normal). The final state leptons in the radiative 
decay mode Bg — s- £+£"7 , apart from having longitudinal polarization, can have two more components of polarization 
( Pt is the component of the polarization lying in the decay plane and is the one that is normal to the decay 
plane) p9| ]. Both Pjv and Pt remain non-trivial for /i^/i^ and t+t^ channel since they are proportional to the 
lepton mass, mi . The different components of the polarization i.e. Pt , Pn , Pt involve different combinations of 
Wilson coefficients and hence contain independent information. For this reason confronting the polarization results 
with experiments are important investigations of the structure of SM and for establishing new physics beyond it. The 
radiative process Bg £^£~"f has been extensively studied in 2HDM and SUSY by various people |^ and the 
importance of the neutral Higgs bosons (NHBs) has been emphasized in the decay mode with /i and r pairs in the 
final state. In this work we study various polarization asymmetries associated with final state lepton (considering 
lepton to be either muon or tau) with special focus on the NHB effects . 
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Bs £+f~7 decay is induced by the pure leptonic decay Bg £^£~ which suffers from hehcity suppression for 
hght leptons {£ = e, n). But in radiative mode {Bg t^£~^ ) this hehcity suppression is overcome because the lepton 
pair by itself does not carry the available four momentum . For this reason, one can expect Bg £~^i~"f to have a 
relatively large branching ratio compared to non-radiative mode despite an extra factor of a. In MSSM, the situation 
for pure dileptonic modes {Bg £^£^ ) becomes different specially ii £ — /i,r and tan/S is large [|[ 01- This is 
because in MSSM the scalar and pseudoscalar Higgs coupling to the leptons is proportional to mftanP and thus can 
be large for £ = and for large tanj3 . The effect of NHBs has been studied in great detail in various leptonic 
decay modes |, || ||, |, |, |, |, [ill ^ . The effect of NHBs on radiative mode Bg £+£--/ has also been studied 
in 2HDM and SUSY [Q. Here we will focus on the NHB effects on various polarization asymmetries within the 
framework of MSSM. 

This paper is organized as follows : In section 2, we first present the Leading Order (LO) QCD corrected effective 
Hamiltonian for the quark level process b — > s£^£~'-f including NHB effects leading to the corresponding matrix 
element and dileptonic invariant mass distribution. In section 3, all the three polarization asymmetries associated 
with the final state lepton are calculated. Section 4 contains discussion of the numerical analysis of the polarization 
asymmetries and their dependence on various parameters of the theory, focusing again mainly on NHB effects in the 
large tan/3 regime. 



II. DILEPTON INVARIANT MASS DISTRIBUTION 

The exclusive decay Bg £+£"7 can be obtained from the inclusive decay b s£+£~7 and further from b s£~^£~ 
. To do this photon has to be attached to any charged internal or external line in the Feynman diagrams for b s£~^£~ 
. As pointed out by Eilam et. al. fl^ , contributions coming from attachment of photon to any charged internal 
line will be suppressed by a factor of /M^ in Wilson coefficient and hence can be safely neglected . So we only 
consider the cases when the photon is hooked to initial quark lines and final lepton lines. To start off, the effective 
Hamiltonian relevant for b s£+£- is |, |||, |, 0, |, | : 

■Heff = ^^^tbVZ I - 2 Jeff ^ sia^.p^l + 75)6 £r£ + ^eff 57^(1 - 75)^ 

+ Cio S7p(l - I5)b + Cq, s{1 + 75)6 ££ + Cq, s{1 + 75)6 } (2.1) 

where p — pi + p2 is the sum of momenta of £^ and £^ and Vtb,Vtg are CKM factors. The Wilson coefficients 
7e//, 9e// and Cio are given in |]l^, |2^. Wilson coefficients Cq^ and Cq^ are given in ^ ^, 0. In addition to 
the short distance corrections included in the Wilson coefficients, there are some long distance effects also, associated 
with real cc resonances in the intermediate states. This is taken into account by using the prescription given in pl| , 
namely by using the Breit-Wigner form of resonances that add on to Qe// : 



there are six known resonances in cc system that can contribute *. The phenomenological factor Ky is taken as 2.3 
in numerical calculations (l^, |2l) . 



Using eq(2.1) we calculate the matrix elements for the decay mode Bg £^£ 7 . When the photon is hooked to 



the initial quark lines, the corresponding matrix element can be written as : 

Ml = ^^VtbV;g { [A e^a^.e* Vg" + ^B [e^pq) - {e*p)q,.)] hH 



'27r ' ' """" ' ' ■ ^ 

+ \C £pa/3.e*Vg" + {e;{pq) - {e*p)q,)] £^15^ } (2.3) 



* all these six resonances will contribute to the channel Bs — > A'^M 7 whereas in the mode Bs — » t+t 7 all but the lowest one J/'I'(3097) 
will contribute because mass of this resonance is less than the invariant mass of the lepton pair (4m|) 
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where A, B, C and D are related to the form factor definition and are define in appendix eqns.(Bl - B4). Here 
and are the polarization vector and four momentum of the photon respectively, p is the momentum transfer to the 
lepton pair i.e. the sum of momenta of and £~ . We can very easily see from the structure of eq.(2.3) that neutral 
scalars don't contribute to A^i . This is due to eq.(B4) given in appendix. 

When the photon is radiated from either of the lepton lines we get the contribution due to Cio along with scalar 
and pseudoscalar interactions i.e. Cq-^ and Cq^. Using eqns. ( p^ - B8) of the appendix |2[ ^ the corresponding 
matrix element is : 



/27r 



(C: 



10 



2mimb 



2mimi, 



2meX 



1 



2piq 2p2<? 



2p2q 2piq 



75^ 



2piq 



(2.4) 



where Pb^ and /b^ are the four momentum and decay constant of meson and pi and P2 are the four momenta of 
£~ and respectively. 

The final matrix elment of Bs t^(r^ decay thus is : 



M = Ml + M2 
From this matrix element we can get the square of the matrix element as : 

^ \Mx\^ + \M2\^ + 2Re{MiM*2) 

with 



(2.5) 



(2.6) 



\M2\^ 



\B 



[p'iipiqf + {P2qf) + 2ml{pqf 



\D 



[p'{(piqf + {p2qf )-2m]{pqf] + 2 Re{B* C + A* D) p\{p2qf - {p^qf)) 



= 4 



V2^ 
1 

iviq) 
1 



(6p2 + 4(p2g)) 



(Ci, 



2mimi, 



Cq,)< 8 + 



1 



(— 2m^_^m^ — niB^p' 
1 



2 .2 , p4 



(2.7) 

2p\p2q)) 



ipiq){p2q) 



+ ( 



1 



(-2m|j to| - m%y + / + 2p^{piq)) + 



(^29) 



(6p2+4(pig)) 



2mirnb 
1 



CqJ 8 + 



ipiqf 



.{QniB^mg + Bnig — tub^P ~ Sm^p +p ~ 8mg{p2q) + 2p {p2q)) 



-i-A0mj + 6p^ + A{p2q)) + 



(6TO3 mf + 8mj — m% p^ — Smfp^ + p"^ — 8nij{piq) 



2p\piq)) 



{P2q) 



(-40m^ + 6/ +4(p2g)) + 



{piq){P2q) 



{Am% ml + lQm\ - IGm^p^ + 2/) 



(2., 



27r 



2memb " { {piq){P2q) 



^^(^) (M)-bl.-p2.) | ^ {ReiB)-^iHpqf-2ipiP2)ipiqf 

{Piqj[P2q) J 2memt [ (Piq){P2q) 



h'2ipip2)ip2qf + 4m^Mq)iP2q)) + Re{G)^P^^^^^^^^ \ (2.9) 

{Piq)[p2q) J 

The differential decay rate of Bs £^-^"7 as a function of invariant mass of dileptons is given by : 



ds ' 2V2^ 16(2^)3 



VtbV,T TTH^ (1 - J) A 1 - ^ A (2.10) 
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with A defined as 
4 



A = § (1 - s)' [ + \B\') i2mj + s) + + \D\'){-Amj + s) 



2 (,'-'10 + 7; l/qJ . 



(1-5)2^1-^ 

64/2 m2 ^1 [ (-1 + I2m2_i6m4_j2)i^(^) _^ (-2s - 8mp + 4s2)a/i" 



32 /s ™| (Cio + tP^CqJ Re{A) ^""^^^ 
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32 (^I^^c-qJ i?e(i?) ^— (2.11) 

2mimb . _ 4m| 



[ (1 - Amj + s)ln{z) ~ 2sd 1 - ^ ] 



where s — p /nin , mj — mj/mg , z = — ^ , are dimensionless quantities 



III. LEPTON POLARIZATION ASYMMETRIES 

We now compute the lepton polarization asymmetries from the four Fermi interaction defined in the matrix element 



eqn.(2.3) and eqn.( p.4| ). For this we need to calculate the polarized rates corresponding to different lepton polarizations. 
These rates are obtained by introducing spin projection operators defined by = ^(1 + JsjS^), where index 
X — L, N, r and corresponds to longitudinal, normal and transverse polarization states respectively. The orthogonal 
unit vectors, Sx, defined in the rest frame of £~' read |19| : 



PL 
Pil 

q X Pi 



V iq X Pi 

S!^ = (0,eT) - (0,e^xeL) (3.1) 
where pi and q are the three momenta of £^ and photon in the center-of-mass (CM) frame of + system. Furthur- 



more, it is quite obvious to note that, S^-Pi — 0. Now boosting all the three vectors given in eqn.(3.1) to the dilepton 
rest frame , only the longitudinal vector will get boosted while the other two (normal and transverse) will remain the 
same. The longitudinal vector after boost becomes ^ : 

Si = f^,^) (3.2) 

We can now calculate the polarization asymmetries by using the spin projectors for £~ as ^(1 + 75 /S*). The lepton 
polarization asymmetries are defined as : 

rfr(g^) _ dr(-s^) 



t this particular choice of polarization is called helicity 
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where the index a; is L, T or A'', representing respectively the longitudinal asymmetry, the asymmetry in the decay 
plane and the normal component to the decay plane. From the definition of the lepton polarization we can see that 
Pl and Pt are P-odd, T-even and CP-even observable while P^ is P-even, T-odd and hence CP-odd observable * 
Our results for the polarization asymmetries are 



Pl(s) 



^ml^ Re{A*C + B*D) ^ 1 - ^ s{l - s)^ 



128/1 m2 



2memh 



CQ2) 



1 



+ {2mj - 8mj ~ s + Smjs - 8mjs + 2rh1f - s^) ln(z) 

1 



4m? 



+ 32 feM (C; 
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^3 /~l \ 

2mfmb (4m^ — s) 



ReiB)\{~s + 3s^)\ll~^^ + 2{mj + mjs ~ s^)ln{z) 



-i?e(C)(l-s)(,sVl 
Re{A) 



irhj 



+ (2m2 - ,s)ln(i) + 32/s,(-^CqJ 



m2(i_s) 



s \ 1 



2mi ln(z) I + Re{D) \s\jl- ^ + (2m^ - s)ln(l) 



/ A (3.4) 



Pt{s) = TT 



2m\^ Re{A*B) - sf - — ^^m^ (Cm 



7^ '^Qi)\T, "-"QJ— q ^ 

Zmirrib Zmiirib (1 — s) 



2mimb 



CqJ i?e(S) 



+ 8/i3. (-^CqJ i?e(A) 



ii-iM±M) + ^e(C) (-2™. + Vi) (1 + 4 
(2mf + Vs) J 



(4m| + s - 12r7i;l + s^) 
i2mi + 



- Re{D) {2mi - VS) {1 - s) 



/A (3.5) 



PAr(s) = TT m£ 



Im{A*D + B*C) (1 - sf J s - irh^ 



8 fs, {Cio + 



-CQ2) 



2mimb ^ (2fhi + ^/s) 



{Im{A){l + s) + Im{D){l~s)} 



+ 8 /s, [tt^CqJ 



y/s- Ami 



2mimb ' {2rhi + \/i) 



{/m(B)(l - s) + /m(C)(l - 8m^ + s)] 



/A 



(3.6) 



with A as defined in eqn.(3.1) and ~ me/niB^ 



IV. NUMERICAL RESULTS AND DISCUSSION 

We have performed the numerical analysis of various polarization asymmetries whose analytical expressions are 



given in eqns.(3.4 - |3j 

Although MSSM is the simplest (and the one having the least number of parameters) SUSY model, it still has 
a very large number of parameters making it rather difhcult to do any meaningfuU phenomenology in such a large 
parameter space. Many choices are available to reduce such large number of parameters. The most favorite among 
them is the Supergravity (SUGRA) model. In this model, universality of all the masses and couplings is assumed at 



-t because time reversal operation changes the signs of momentum and spin, and parity transformation changes only the sign of momentum 
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the GUT scale. The minimal SUGRA (mSUGRA) model has only five parameters (in addition to SM parameters) 
to deal with. They are : m (the unified mass of all the scalars), M (unified mass of all the gauginos), tan(3 (ratio 
of vacuum expectation values of the two Higgs doublets), A (the universal trilinear coupling constant) and finally, 
sgn{fi) § . 

It has been well emphasized in many works |^, |l^ that it is not necessary to have a common mass for all the 
scalars at the GUT scale. To have required suppression in K'^ — mixing, it is sufficient to have common masses of 
all the squarks at the GUT scale. So the condition of universality of all scalar masses at the GUT is not a very strict 
one in SUGRA. Thus we also explore a more relaxed kind of mSUGRA model where the condition of universality of 
all the scalar masses at the GUT scale is relaxed with the assumption that universal squark and Higgs masses are 
different. For the Higgs sector we take the pseudo-scalar Higgs mass (m^) to be a parameter. Over the whole MSSM 
parameter space we have imposed a 95 % CL bound , consistant with CLEG and ALEPH results : 

2 X 10""* < Br{B ^ Xs-f) < 4.5 x 10"^ 

Figure(|l|) shows plots of the differential Braching ratios of Bg £^£^j for Icptons to be /i and r . The prediction 
of the Branching ratios for Bg £^£^'y are : 



TABLE I: Branching ratios for Bs ~* i+l'j {£ = ^, r) 



Model 


Br(B, ^ ^i'^^^^"f ) 


Br(B, ^ 




Standard Model 


5.53 X lO"** 


6.57 X 


io-« 


mSUGRA " 


6.86 X 10"* 


3.59 X 


10-^ 


SUGRA " 


1.21 X 10"^ 


1.31 X 


10"® 



"The mSUGRA and SUGRA parameters are defined in Figure{y). These values are of the same order as estimated by Xiong et. al. ^ 

We have plotted various polarization asymmetries {Pl , Pt and Pjv ) in the three models - SM, mSUGRA and 
SUGRA in Figures (^,||,||) for Bg ji^ and Bg t^t^^ as a function of s (scaled invariant mass of the dilepton 
pair). 

Now we try to analyse the behavior of the polarization asymmetries on the parameters of the models chosen 
(mSUGRA, and SUGRA). For this analysis we consider the polarization asymmetries at dilepton invariant mass (s) 
away from the resonances (the J /^) resonances (we choose s = 0.68 for our analysis) . The main focus of the analysis 
is NHB effects on polarization asymmetries. These effects crucially depend on tanP and pseudoscalar Higgs mass 
(toa). 

In mSUGRA model the Higgs mass (at electroweak scale) depends crucially on the universal mass of the scalars 
and tanj3 . To illustrate this crucial behaviour, we have plotted various polarization asymmetries as a function of 
tan/3 for different values of unified scalar mass (m) in Figs. (||, H) . As can be seen from these figures, Pl shows 
large deviations from the SM values and over a significant portion of the allowed region, even shows a sign flip, 
provided ton/3 is sufficiently large. Similar behaviour is also there for Pt ■ On the other hand, the predictions for Pjv 
don't differ substantially from SM results but the mSUGRA predictions can change Pn by more than 50 % with an 
appreciable increase in tan/3 . 

For SUGRA model we have plotted ( Figs. (^, |^, |l^) the polarization asymmetries as a function of pseudoscalar 
Higgs mass (m^) for various values of tan/3 . In SUGRA we expect more variation of all the polarization asymmetries as 
compared to their SM values because here we have Higgs mass (pseudo-scalar Higgs mass) as an additional parameter 
along with tan/3 . As we can see from Figure(^) the variation of Pl is more substantial in SUGRA model. In fact 
for fairly large region of SUGRA parameter space, Pl can be opposite in sign as compared to SM case . Pt can vary 



S our convention of the sgn{^) is that fj. enters the chargino mass matrix with positive sign 
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upto five in magnitude when compared with the SM value over the large region of allowed parameter space and for 
the parameter space we have taken into consideration, the predicted value of Pt in SUGRA is opposite in sign to the 
SM value. Again does not show as much deviation as observed for Pj, and Pt but the variation can still be upto 
an order in certain region. 

Summarizing the results of the numerical analysis : 

1 . From Figures (^|j|) it is clear that the longitudinal and transverse polarization asymmetries (P^ , Pt ) can have 
substantial deviation from their respective Standard Model values over the whole region of dilepton invariant 
mass (s), while Figure(^) indicates deviation for P/v from SM values for a limited region of the dilepton invariant 
mass. 

2. As we have pointed out earlier |^ that for the inclusive process B — s- Xs£~^£~ there is not much deviation from 
SM results in mSUGRA model. But in radiative dileptonic decay mode, mSUGRA predictions also show large 
deviations (at least of Pl and Pt ) from SM results, making it possible to use polarization asymmetries to test 
mSUGRA model. This is mainly because in the bremmstrahlung part of the matrix elment {M2), the Wilson 
coefficient Cq^ adds on to Cio via the combination (Cio + 2mfmb ) which effectively increases the SM value 
of CiQ. This doesn't happen for the process B Xst^i~ and this numerically is the reason for the scalar 
exchanges affecting the Bg £+^"7 process more than the semi-leptonic one. 

3. From Figs. (^, ^, ^) we can see that the polarization asymmetries show a general enhancement with increase in 
tan(3 and they decrease as the universal scalar mass (m) is increased. This is expected because the Higgs boson 
mass increases with ni and thus the contributions of scalar (Cq^) and pseudoscalar (Cq^) type interactions 
decrease. 



4. As can be seen from the structure of the analytical expressions for various polarization asymmetries ( eqn.(B.4 



3.5, |3.6| ), they are all different analytic functions of various Wilson coefficients and hence contain indepen- 
dent information. These asymmetries, hence, can also be used for accurate determination of various Wilson 
coefficients. 

In conclusion, we can say that the observation of the polarization asymmetries can be a very useful probe for finding 
out the new physics effects and testing the structure of the effective Hamiltonian. 
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APPENDIX A: INPUT PARAMETERS 



ms, = 5.26 GeV , nic — 
■m^ = 0.106 GeV , rrir — 
mw = 80.4 GeV , rrit = 
Gf = 1.17 X 10-5 GeV"^ , a = 



1.4 GeV , m, = 0.2 GeV 
1.77 GeV , TOf, = 4.8 GeV 
176 GeV , \VtbVt%\ = 0.045 
: jig , T{mB^) — 1.6 X 10"-^^ sec 



APPENDIX B: 

Definition of A, B, C and D defined in eq(^.3|) are : 

A = ^ mfG,{p') - 2Te//^G2(/)], 



B 



D = 



C'lo 



Flip'). 



where the form factors definition chosen is 



(7| S7m(1 ± 75)6 \Bs) - -r- {e,o.f3.elppq,G,ip^) T ^[{e;{pq) ~ (e*p)gJ]Fi(p2)} 
(7I siaf,^p^{l ± 75)6 \Bs) = — ^ {ef,ai3ae*aPi3qaG2ip'') ± i[(e;i(w) - (e*P)9p)]-p2(/)} 



multiplying eq. ( B2 ) with p^ and using equation of motion we can get relation : 

(7I s(l±75)& \Bs) - 
The dcfination of form factors we are using for numerical analysis is p2[ ; 



Clip' 



1 



GeV , G2{p' 



3.74 



1 -PV5.62 
Fi(p') = , Gey , F2{p') 



1-392/40.5 
0.68 



GeV 



GeV 



1-P76.52 ' ' l-p2/30 

Identities used in calculation of matrix element when photon is radiated from lepton leg : 

(0| sb \Bs) = 
(0| S(T^.(l + 75)fo \Bs) = 

(0| 57^756 \B,) - - i/s.Ps.p 
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FIG. 1: Branching ratios for Bs —^ t^i 7 with £ = fj, (above) and £ = t (below). mSUGRA parameters are : m = 200 GeV, 
M = 450 GeV , A = , tan /? = 40. Additinal parameter for SUGRA (the pseudo-scalar Higgs mass) is taken to be ruA = 306 
GeV 
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FIG. 2: Pl for Bs -» £+£"7 with e = fx (above) and £ = t (below). mSUGRA parameters are : m = 200 GeV, M = 450 
GeV , A = , tan /3 = 40. Additinal parameter for SUGRA (the pseudo-scalar Higgs mass) is taken to be niA = 306 GeV 
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FIG. 3: Pl vs ton/3 at J = 0.68 for -> r+r 7 in mSUGRA model, other parameters are : M = 450 GeV , ^ = 0. 
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FIG. 4: Pl vs jtia at s = 0.68 for t+t~^ in SUGRA, other parametes are : m = 200 GeV , M = 450 GeV , ^ = 0. 
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FIG. 5: Pn for Bs ->■ i+Tj with £ = fj, (above) and ^ = r (below). mSUGRA parameters are : m = 200 GeV, M = 450 
GeV , A = , tan /3 = 40. Additinal parameter for SUGRA (the pseudo-scalar Higgs mass) is taken to be niA = 306 GeV 
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FIG. 6: Pat vs tanP for t+t'j at s = 0.68 in mSUGRA , other parameters are : M = 450 GeV , A = 0. 
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FIG. 7: Pat vs mA at s = 0.68 for Bs -> T+r"7 in SUGRA , other parametes are : m = 200 GeV , M = 450 GeV , A = 0. 
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FIG. 8: Pt for Bs t+T^ with I = ^, (above) and £ = r (below). mSUGRA parameters are : m = 200 GeV, M = 450 
GeV , = , tan /3 = 40. Additinal parameter for SUGRA (the pseudo-scalar Higgs mass) is taken to be ruA = 306 GeV 
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FIG. 9: Pt vs ton/3 for Bs -> T+T~-y at s = 0.68 in mSUGRA , other parameters are : M = 450 GeV , A = 0. 
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FIG. 10: Pt vs tua for S„ t+t~j at s = 0.68 in SUGRA , other parameters are : m = 200 GeV , M = 450 GeV , A = 0. 



